Conical soliton escape into a third dimension of a surface vortex 
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We present an exact three-dimensional solitonic solution to a sine-Gordon-type Euler-Lagrange 
equation, that describes a configuration of a three-dimensional vector field h constrained to a sur- 
face p-vortex, with a prescribed polar tilt angle on a planar substrate and escaping into the third 
dimension in the bulk. The solution is relevant to characterization of a schlieren texture in nematic 
liquid-crystal films with tangential (in-plane) substrate alignment. The solution is identical to a 
section of a point defect discovered many years ago by Saupe [Mol. Cryst. Liq. Cryst. 21, 211 
(1973)], when latter is restricted to a surface. 

PACS numbers: 



I. INTRODUCTION 

Topological defects are central to a complete descrip- 
tion of ordered phases of condensed matter, ranging from 
superconductors to liquid crystals^. Defects' energet- 
ics controls the stability of the ordered state to ther- 
mal fluctuationaS.!^, random material heterogeneities^ 
and external perturbations^. A complete rigorous 
classification^ is now available for most hulk ordered 
states. 

This has been particularly fruitful in understanding a 
rich variety of topological defects that are found in liquid- 
crystal phases. However, in many physical contexts, as, 
for example, arising in liquid crystals confined inside a 
thin display cell, much of the physics is controlled by 
a substrate interaction which competes with the bulk 
energetics^. In such surface-dominated situations, only 
an incomplete understanding of defects structure and sta- 
bility is available. 

One important and extensively studied example of this 
type discovered by Meyer— is that of a uniaxial nematic 
liquid crystal confined to a thin long capillary with a 
homeotropic alignment at the cylindrical surface. The 
resulting boundary condition forces an integer winding 
of the nematic director field, which for a two-dimensional 
{xy) field would trap a vortex line along the axis of the 
capillary. However, such defect is unstable for a three- 
dimensional (3D) director field and away from the bound- 
ary exhibits an escape into the third dimension, removing 
the line singularity as described by Meyer's solution^. 

A familiar schlieren surface texture seen in phase con- 
trast microscopy is a hallmark of nematic liquid crys- 
tals, reflecting surface-induced disclinations (vortices in 
the nematic director field). The texture details, e.g., ap- 
pearance of integer versus half-integer vortices have been 
suggested to distinguish between the uniaxial and biax- 
ial nematic states^i^. New advanced bulk imaging tech- 
niques, such as, for example, the fluorescence confocal 
polarizing microscopy, have also allowed imaging of the 
full three-dimensional textures associated with such sur- 
face defects^". 

Motivated by the above discussion, here we consider a 
problem of a 3D nematic with a planar (parallel) surface 
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FIG. 1: A nematic director field n(x) constrained to a 27r- 
vortex on a surface z = Q with a planar alignment rendered 
in 3D in (a) and in a 2D projection in (b) . The corresponding 
bulk texture configuration n(x) that minimizes the Frank en- 
ergy is calculated analytically in this paper and is illustrated 
in Fig.©. 



alignment, with an integer vortex imposed on a substrate, 
as illustrated in Fig.ll]). In contrast to the long capil- 
lary case^., that clearly exhibits translational invariance 
along its axis, reducing it to one dimension (ID), here the 
system is manifestly three-dimensional [two-dimensional 
(2D), once azimuthal symmetry is included], and there- 
fore in principle considerably more complicated. 

Here we present a derivation of an exact solution to the 
single Frank elastic constant Euler-Lagrange (E-L) equa- 
tion that describes a bulk texture induced by a surface 
27rp-vortex, with p the integer azimuthal vortex wind- 
ing number. It is described in terms of the polar angle 
9{r^, z) = 9s{z/r±) of the director field n(x) that we find 



to be given by 

9s{t) — 2arccot 



(t+vt^ny 



(1) 



Illustrated in Fig. ([2]) for p = 1, the bulk texture is a 
conical soliton giving the nematic director's escape into 
the third dimension away from the imposed surface 2'Kp- 
vortex with a strong planar alignment. 
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FIG. 2: Exact conical soliton solution (here projected onto 
z ~ rx_ plane) of the nematic director field, n(x) , describing 
escape into the third dimension of a surface 27r-vortex. 

The director n(x) configuration corresponding to 
Eq.(IT]) coincides with a section of the Saupe's point -p- 
defect*^ [tan 6^/2 = (tan(5/2)IPl, with (5 as a polar angle 
of the spherical coordinates] when restricted to a sub- 
space above a planar substrate. For p = 1 the above 
result reduces to a well-known simple texture, given by 
half of the hedgehog (skyrmion) configuration, n(x) = x. 
Application of Saupe's p = 1-defect to a surface vortex 
problem was also previously explored by Klemar>i^ and 
was shown to satisfy the simplest homogeneous boundary 
conditions arising from a model pinning potential^^. 

To summarize our contributions, we present a deriva- 
tion (in cylindrical coordinates, mapping the Euler- 
Lagrange equation to that of a dissipative particle with a 
time-dependent mass) of Saupe's solution [Eq.([T])]. This 
approach is likely extendable (even if approximately) to 
a study of other interesting problems of surface defects, 
where the reduction to Saupe's point defect solution no 
longer holds. Our slight generalization of Saupe's solu- 
tion allows us to discuss and connect to weak and strong 
anchorings. Finally, we present an analysis of the en- 
ergetics, comparing the conical soliton escape to other 
competing textures. We discover a counterintuitive de- 
pendence of the p-vortex energy on p, showing that it 
is asymptotically linear in p, in contrast to the standard 
]? dependence. This has important implications for the 
stability of p > 1 surface vortices over their fission into p 
lower winding (p = 1) vortices. 



II. MODEL 

We consider a 3D model of a nematic liquid crystal, 
characterized by a nematic unit director field, n(x), with 
X = (rj^, z). The energy is given by a Hamiltonian 

H^H,i + Hs, (2) 

where H^i is the bulk elastic energy of the Frank model 

1 



H, 



d V I dz 



Ki{\I ■ hy + K2[n ■ {\I X n)Y 



-|-i^3[n X (V X n)f 



(3) 



and Hs is the surface pinning energy, localized at z = 0, 
Hs = j d\^dzVs{vi_)5{z){z ■ hf, (4) 



with z as the surface normal. In the single elastic con- 
stant approximation, Ki = K2 = K^ = K, the elastic 
energy reduces to 



K f 



(5) 



In above we have dropped the boundary terms as they 
do not affect the E-L equation. Parametrizing the unit 
director field 



n = (sin 6 cos 0, sin 9 sin 0, cos ( 



(6) 



in terms of polar and azimuthal angles 9 and (j), H^i re- 
duces to 



Hei = Y I d^^^dz 



{W9y + sm'^ 9 {W(l>y' 



(7) 



We can include surface pinning through a boundary con- 
dition on n(r±, z = 0) = fio(rx), finding the correspond- 
ing solution and then minimizing over no{r±) in the pres- 
ence of T4(r±). 

We focus on the solution 7i(x) subject to a constraint of 
a 27rp-vortex {p £ Z) at z — 0. The 2Trp surface winding 
is imposed by taking 



(rj.,z ^ 0) ^ pip + ipo, 



(8) 



where ip — arctan(j//x) is the azimuthal angle of the 
cylindrical coordinate system x — {r j_ cos ip,rj_ sirup, z). 
The arbitrary constant angle (po gives a family of textures 
induced by spiral surface defects for < (po < -k 12. These 
interpolate between a pure splay surface "aster" defect 
for ipQ = 0, [illustrated in Fig.([T])(a)] and a pure bend 
surface "vortex" defect for 930 — 1^/2. 

The resulting elastic energy is then given by 



K I , 



{V9f + ^ sin^ 



(9) 



leading to the Euler-Lagrange equation that determines 
the texture configuration ^(x), 



V^6'--^ sin 261 = 0. 
2r2 



(10) 



Focusing for siniphcity on azimuthally symmetric 
boundary conditions, we search for a (^-independent so- 
lution 9{r,z), satisfying 

r^dle + rdrO + r^dle - ^ sin 26* = 0, (11) 

with a surface constraint 9{r,z = 0) = do{r). We have 
simplified the notation by denoting r± = r. 



III. SOLUTION OF THE EULER-LAGRANGE 
EQUATION 

Despite the fact that the E-L equation, Eq. pT|) is non- 
linear and two-dimensional, its one- (and periodic array-) 
soliton solution can be found cxactly^^. The intuition for 
the form of the solution can be obtained by neglecting 
the r-derivative terms and then noting that the resulting 
equation is of a standard ID sine-Gordon type along z, 
with a period tt. It thus admits a soliton solution con- 
necting tilt angle 9{r, z = — oo) = tt to Q{t, z = -|-oo) = 
with the soliton width at the transverse distance r from 
the vortex given by i^zir) — """■ 



A. Exact conical soliton solution 

Motivated by the above observation and by the trans- 
lational invariance of the E-L equation along z, we search 
for a soliton solution of the form 



[r,z) 



z + Zq 



(12) 



We note that this restricted form precludes a study of 
other than a constant boundary condition at z = — zq. 
Since in general the symmetry dictates a nontrivial ra- 
dial variation in the director tilt angle at the surface with 
the distance r from the vortex, we anticipate that the 
above form of the solution is an exact description only 
for an infinitely strong planar alignment on a substrate 
aX z — —zq = 0. For a finite planar surface anchoring, 
given large azimuthal strain near the vortex, we expect 
a meron configuration with 0o(^) ~ in the vicinity of 
the vortex (near r = 0) and growing to 7r/2 with increas- 
ing distance r from it. As we will see below, the radial 
surface variation 0o(^) = S{r, z = 0) can be qualitatively 
captured by the solution 9s{zo/r) at z = 0, by adjust- 
ing zp. While this single degree of freedom (zq) is in 
principle insufficient to capture an arbitrary form of the 
surface boundary condition, OqIt), we proceed to explore 
this class of solutionsi^. We expect it to be a good ap- 
proximation for strong planar anchoring, characterized 



by a vanishing zq and corresponding to 9o{r) « tt/2 for 
nearly all r, excluding a small core region of radius zq. 
We will treat the weakly anchored case in a complemen- 
tary way. 

As a function of the scaling variable t = {z + ZQ)/r, the 
E-L equation simplifies to 



mit)9, +l{t)e,-^ sin 29 s = 0, (13) 



where 



m(i) = i^ + 1, 
7(i) = t. 



(14) 
(15) 



Solutions of Eq. p^ can be most easily obtained by 
its identification with Newton's equation for a particle 
at position 9{t) at time t, moving in a periodic potential 

2 

V{9) = ^ cos 29 and characterized by time-dependent 
mass and friction coefficients, m(t) and 7(t), respectively. 
This type of identification is quite analogous to a stan- 
dard sine-Gordon model, where, in contrast, the fictitious 
particle has a constant mass and no friction. In this lat- 
ter case the solution is easily obtained by a guaranteed 
existence of an integral of motion, energy of the parti- 
cle, which reduces the solution to a single integral. In 
our problem the time dependence of the mass and finite 
friction at first sight would be expected to preclude the 
existence of such "conservation of energy" integral of mo- 
tion. However, a key observation is that the two effects 
can exactly compensate each other if the condition 



1 



lit) = -m{t) 



(16) 



is satisfied (as it is in our problem) , and leads to an "en- 
ergy" conservation law 



d_ 
dt 



-m(i)0'/ + ^cos20. 



0. 



(17) 



Indeed this is, guaranteed by the fact that the E-L 
Eq. ljlOp came from a minimization (of Hpi) principle. 
The resulting integral of "motion," 



1 2 2 

-m{t)9'l + — cos 261, = —E, 

2 4 4 



(18) 



with p^E/A as the fictitious particle's energy, then easily 
allows us to reduce the problem to a single integral 



to \/m{t') 



= -V2 



d9' 



V2 V^- cos 26*' 



(19) 



where in above, we have made a choice of the negative 
square-root. The parameter to defined by 9s{to) = 7r/2 
determines the tilt angle 9q at the boundary at z = — zq. 
The other constant of "motion," E, is also crucial to the 
nature of the solution. It is quite clear that for a half- 
infinite space boundary conditions (see Fig. ([I])) E must 
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FIG. 3; Solitonic p = 1 solution from Eq.((2T]) for (a) to = 0, 
corresponding to a perfect planar alignment, 6o = 7r/2 at 
z = — zo, and (b) to = —0.5, corresponding to a uniform tilt 



of ( 



1.11 radians at the z = —zq boundary. 



be chosen to be i? = 1~ so that the solution is a single 
soliton in t. In the mechanics analogy it corresponds to 
a particle at t = — oo, starting out at = tt, with the 
energy just equal to the potential energy, i.e., with an 
infinitesimally vanishing initial velocity, rolling down the 
hill during — oo < t < to and then climbing back up to 
the top at 6* = as t ^ +oo. 

Using ■m{t) = t^ + 1, Eq.([T9]) is easily integrated, 



P 



pin 



dt' , 

t + VWTT 



<o + V^§ + i 



dO' 
/^/2 sin 6*'' 

= lncot(6's/2). 



and leads to our main result 



9s{t) = 2arccot 



t + Vt2+T 



h + ^M+l 



(20) 



(21) 



illustrated in Fig. ([3]). The corresponding director field 
texture is illustrated in Fig. ([J). The solution satisfies the 
Euler-Lagrange equation with a uniform tilt 6*0 boundary 
condition. 



9{r, z = —zq) = 9q = 2 arctan 



to + 



\/^ 



(22) 



FIG. 4: Surface tilt angle 6o(r) for a (p = 1) 27r-vortex, show- 
ing surface escape into the third dimension (a meron), con- 
fined to a radius zo. 



on the z — ^zo surface. Clearly, ^o is also the asymptotic 
tilt angle at large r (vanishing t), and in terms of it, the 
solution can be equivalently written as 



j(t) — 2arccot 



cot 



|(<+v^ 



+ 1 



(23) 



This is illustrated for p = 1 in Fig. ((31). 

We focus on the asymptotically planar alignment, 
9{r -^ oo^z) — 9o — 7r/2, corresponding to to = 0. On 
the physical surface boundary at z = 0, the tilt angle is 
then given by a nontrivial function of r, 



9ir,z = 0) = 9o{r)=9s{zo/r), 



= 2arccot 



za/r+ yz^/r^ + 



(24) 
(25) 



that describes the escape into the third dimension [van- 
ishing 9o{r)] on the surface z = inside a disk of radius 
zq, as illustrated in Fig. ([3]). As anticipated above, zq 
allows only a single parameter adjustment of the bound- 
ary condition, physically controlled by Vj^. An infinitely 
strong surface anchoring, Kj ^ oo, gives a perfectly pla- 
nar alignment, Oo^r) — 7r/2, characterized by zq — > 0. 

We note, however, that a more general boundary con- 
dition, 9o{r) can be imposed by generalizing above exact 
solution to an r-dependent ZQ{r). Although the resulting 
9s{[z + zo{r)]/r) is no longer an exact solution to the E-L 
equation, for a small drZo it is an accurate approximation 
and can be employed as a good variational ansatz. 



B. Pontryagin index of conical solitons 

The solitonic solutions 9p{x.) ~ 9s{t), Eq. (P5)) (indexed 
by p), together with 0p(x) — pip, give the unit director 
field np(x) according to parametrization Eq.®. When 
restricted to a two-dimensional closed surface, e.g., a 
sphere Sf in coordinate space x, np(x) gives a map- 
ping of this coordinate sphere 5*2 into another sphere S'J, 



where n "lives." Such mappings fall into topologically 
distinct classes, that form the second homotopy group, 
H2iSs) — 2, corresponding to distinct ways of wrapping 
a coordinate sphere around a target space sphere. The 
classes are characterized by the Pontryagin topological 
index 



Q = T— dakei-jkfi ■ {difi x djfi), 
87r J 



(26) 



where dak is the fcth component of the infinitesimal sur- 
face element pointing along the local surface normal. We 
have computed Q for our director field solutions np(x), 
and found that Q = p. 



C. Weak surface pinning (small 9) analysis 

For weak surface pinning the tilt angle 9o{r) is small, 
corresponding to a large zq, and for a large range < 
r < zo, solution Eg. ([23]) reduces to 



6{t) 



Oo 



(f + VFTT) 



i-t 



(27) 



We can compare this result with that obtained by 
studying the linearized^"* Euler-Lagrange equation, 

r^dle{r,z) + rdrd{r,z)Wdle{r,z)~p^e{r,z) = 0. (28) 

This differential equation is separable. Letting 9{r, z) = 
R{r)Z{z), it becomes 



^ R ^ r R r^' Z 



(29) 



where the sign of the constant — fc^ is chosen negative 
to ensure a well-behaved solution that decays at large z. 
Keeping only the decaying solution, Z{z) is given by 

Z{z) = Zoe-'''. (30) 

Noting that the equation for R{r) is the Bessel equation 
of order p in variable kr, the full solution of the E-L 
equation for weak pinning is given by 



9{r,z)~ / dkakJp{kr)e 
Jq 



-kz 



(31) 



where the coefficients ak are determined by the boundary 
condition at z = 0, namely by OQ{r) = d{r, z = 0). Using 
the orthogonality of Bessel functions, these are given by 



ak = k drr9{r,0)Jp{kr). 

Jo 



(32) 



We note that, in contrast to the full solitonic solution, 
Ea. (j2ip [where we were only able to impose a boundary 
condition with a specific r dependence, Ooir) = 0s{zQ/r), 



displayed in Fig. ([4])], here we can obtain a solution d(r, z) 
for an arbitrary r-dependent boundary condition 9(r, z = 

0) - eo{r). 

To compare to the full solution, we choose a constant 
boundary condition Oq, for which 



ak = k drr9oJp{kr), 



= 9op/k, 



(33) 
(34) 



-O+fer 



where a convergence factor e '' had to be introduced 
to make the integral into /„ dxxJp{x)e~^ ^, which is 
well defined and equal to p. Using these expansion coef- 
ficients ak inside Eq. (j3ip we find 



1"°° dk 
r,z) - 9opJ ^Jp{kr)e-''^, (35) 

(36) 




in complete agreement with the small 6*0 limit [Eq. (|27|) ] 
of the full solitonic solution. 



IV. ENERGETICS 

A. Conical soliton energy 

The bulk elastic energy corresponding to the surface 
vortex solution, 9s{t), found above is straightforwardly 
computed by plugging into the elastic Hamiltonian and 
evaluating the spatial integrals. We thereby obtain 



El,''\9o,zo) = H,i[9,{t)], 

= -K I d^rdz 



(V^.)^ + ^sin^ 



(37) 
, (38) 



= 2ttp^K / dr / dt sin^ 9s, (39) 

Jo Jza/r 

where we took advantage of the energy integral of "mo- 
tion," Eq. p^ to eliminate (VSg)^, and Lr is the extent 
of the system in the radial direction. Using the explicit 
solution for 9s{t) and defining 



:{t) = t+Vf^ + l, 
xo = a:(io) = (tan(0o/2))'/^ 



(40) 



we obtain 



L,. 



EiP\9o,zo) = 2ttK dre{xo,x{zo/r),p), (41) 

"'0 



where 
e{xo,x{zo/r),p) = 2p^ 



dx- 



x(zo/r) x2{{x/xo)P + ixo/x)P) 

(42) 



s(eo) 




OJ 



e. 



e(p)/p 




FIG. 5: Energy e{6o, 1, 1) (in units of 2-KKLr) of the soliton 
texture due to a 2-k surface vortex as a function of the surface 
tih angle 9o, together with its parabolic approximation, e ~ 
a^Oi with a fitted to be 0.64 (appearing indistinguishable). 



FIG. 6: Energy e(6'o, l,p)/p [in units of 2TvKLre{eo, 1, 1)] of 
the soliton texture due to a 2np surface vortex as a function 
of its topological charge p, displayed for surface tilt angles 

do = 7r/2, 7r/3, 7r/4, 7r/6 (top to bottom). 



1. p=l vortex energy 

Specializing to the case of p = 1 surface vortex, above 
energy is simplified and can be calculated analytically 



E^^H 



L,. 



h,zo) = -iTiK / dr I dxxl—^ 



1 



x(za/r) 



r2^2■ 



L,. 



2ttK dre{xo,x{zo/r),l), 

Jo 

/■OO 

2t:Kzo / dtt^^e{xo,x{t)), 

Jzo/Lr 



(43) 
(44) 



where 

e(a;o,a;,l) 



Substituting this into Eq. p^ we find 

EP{7r/2,zo) = 27TKLr9{zo/Lr) 
where the scaling function g{zQ) is given by 

5(^0) = 2 + arccot f zo + \J zl + 1 j 

— arctan I zq + y^o + -'^ ) + n^o ^^ 

2. p vortex energy 

Ap) 



(50) 



1 + Zn 



m 



:Xo 



xo arccot 



For p > 1 charge vortex, energy El can only be eval- 

xq 1 Xo uated numerically. Focusing on zq = for simplicity, 
1 arctan — . 

X Xo X 



(45) 



E^/'>{eo,0) = 2TTKLre{xo,l,p), 



(52) 



For a vanishing zo, x — x{zo/r) = 1, and we find 

£;«(0o,O) = 2TiKLre(xo,l,l), 
where 



£(a;o,l,l) = -Xo + ^——i^ - Xo m:ccotxo 



where we evaluated e(xo, l,p) numerically and displayed 
e(a;o, l,p)/p as a function of charge p for various values of 
Xo in Fig. [51 As can be seen from this figure, despite the 
(4o) fact that the naivep dependence of e(a;o, l,p) in Eq.[42lis 
the standard p^ found in a 2D vortex, the p dependence 
coming from the integral reduces it to an asymptotically 
linear one at large p. 



•-0 



— arctanoio 



(47) 



8 2 



+ 2(2:0-1), 



for Xo < 1 

for Xo -^ V 



Using the relation 9o[xo) [Ea. (|^ ] to express the soliton 
energy in terms of the surface tilt angle 9o, we obtain 
the energy Eg (6*0,0) of a surface 27r- vortex plotted in 
Fig.®. 

For asymptotic planar alignment 6o{r -^ 00) = 7r/2 
(to = and xo = 1) 



e(xo, l,p> 1) w (1 - cos6'o)p. (53) 

This finding has important qualitative implication that 
(48) °^^ winding-p surface vortex (one p-boojum) has a lower 
energy than p winding-1 surface vortices {p 1-boojums). 
This contrasts strongly with the standard 2D vortex case 
where a p- vortex always has a higher energy than p unit- 
vortices and thus always fissions into them. 



e(l,x) 



TT 

2 



arctanx + arccotx. 



(49) 



B. Competing states 

We can compare the energy of the solitonic state, 
9s[{z + zo)/r\, discussed above with competing states il- 



(a) 



(b) 
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FIG. 7: Competing states for a surface 27r-vortex boundary 
condition, with (a) vortex line and (b) domain-wall texture 
extensions into the bulk. 



lustrated in Fig. ([7]). To this end, we estimate energetics 
by simple scaling analysis for a system of size Lr x L^, 
focusing on the strong planar alignment, 6q = 7r/2. 



1. Vortex line 

One competing configuration is the "vortex line" state 
that extends the surface 27r-vortex into a straight vortex 
line with 0{z) = 9o = 7r/2, independent of z. The energy 
of such a state is clearly given by 



E., 



vortexline 



K 

y 



-^2'Krdrdz oc KL^ In — , (54) 



where a is the core radius of the vortex line, set by the 
coherence length. 



2. Domain wall 

Another possible texture is that of a 2D domain wall 
of thickness a, where 9{z) exhibits a uniform (i.e., r inde- 
pendent) escape into the third dimension, changing from 
^0 = T^/'^ to = over a microscopic distance a. The 
corresponding energy is given by 

K f^- , r , /7r/2 

Edomainwall = TT / ^TTrdr / dz 

2 Jo Jo \ a 

scaling with the area of the cell. 



ex KLl/a 



(55) 



3. Conical soliton surface vortex 

The energy of a conical soliton surface vortex can be 
similarly estimated. We first note that by virtue of the 



E-L equation, all three (z derivatives, r derivatives, and 
sin^ 0) contributions are comparable, and therefore we 
can focus on one of them. Estimating the elastic energy 
based on the z derivatives, we observe that the strain 
is confined to a soliton width along z that at radius r 
is given by ^2 « r. Thus the estimate is quite similar 
to the previous case of the domain wall but with strain 
spread out over region between the cones z = r and z = 
rather than confined to a slab < z < a. This leads to 
an estimate 

L^ pr / _ /o\ 2 



Es^K 



2TTrdr 



dz 



■k/2 



ex KLr 



(56) 



which agrees qualitatively with our exact computation 
Eq.®. 

Since the conical soliton solution scales only linearly in 
Lr, we conclude that the domain- wall solution (scaling 
as L^), is not competitive with the other two solutions. 
On the other hand, the relative competition between the 
vortex line and conical soliton solution depends on the 
relative ratio of L^ and L^-. 

For Lz > Lr clearly vortex line is energetically more 
costly and conical soliton texture is the preferred state. 
On the other hand for a thin cell with width L^ = w < 
Lr a more detailed analysis is required. The vortex line 
energy is still clearly given by EyortexUne = Kwln ^. 

To compute a conical soliton energy in a cell of a finite 
width w requires an extension of the solution to a finite 
geometry. For a finite width cell with free and planar 
boundary conditions on the top and bottom substrates, 
respectively, our exact solution, 9s{t) is a good descrip- 
tion. Its energy can be simply estimated. Examining 
Fig- ©J it is clear that for w < Lr, there are two addi- 
tive energy contributions of this texture. For the region 
< r < ■«;, the contribution is identical to that made in 
Eg. ([55]) . On the other hand, for region r > w, the strain 
field is that of a 27r-vortex line with length w and core 
radius w. Putting these two contributions together, we 
find 



E, 



conic— soliton 



Kw + Kw\n—, (57) 

w 

« Kw{l + In — ^ - In — ) < Eyortex-line, 

a a 

(58) 

for cell thickness w ^ a. 

Unfortunately, we have been unable to find an ex- 
act solution for the experimentally more relevant case 
of non-free (e.g., symmetric planar) boundary conditions 
on both substrates. The difficulty has to do with the fail- 
ure of a periodic soliton solution [obtained by picking the 
integration constant in Eq. p^ to be _E < 1 and match- 
ing its period to the width of the cell; see the Appendix] 
to enforce fixed z (as opposed to fixed t) boundary con- 
ditions. 

However, a good approximate symmetric solution, il- 
lustrated in Fig.®, is given by 



^ir,z)=9s{{w/2-\z\)/r). 



(59) 
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FIG. 8: A solitonic texture describing escape into a third di- 
mension of two 27r-vortices confined to top and bottom sub- 
strates of a finite width cell. 



Describing a cell with two planar aligning substrates at 
z = ±w/2, its only shortcoming is a small slope discon- 
tinuity in the z derivative at z = (the center plane of 
the cell). 



Thus our solution can be expressed in terms of the Jacobi 
elliptic function sn((/), k) defined by 

sn[F{(t),k),k] = sin0. (61) 

sn(i, k) is an odd periodic function resembling a 
smoothed out square wave. For fc > 1 it interpolates 
between a single soliton for /c = 1+ (half a period of a 
square wave) and k^^sinkt for k ^ 1. For fc > 1 the 
period of sn(t, fc) is given by 2F[sin^^(l/fc), fc]. 

To establish a direct relation we change variables </)' = 
e' - 7r/2, finding 



F(0,fc) 



sn (sin0, fc), 



7r/2 ^/T 

fc L/2 



k? cos^ 4)' 

dO' 

^/Ek - cos 26' 



(62) 
(63) 

(64) 



where Sfe = (2 - k^)/k'^. In this notation, Eq.([Tn]) be- 
comes 



dt'^= 

to Vt^ 



-V2 



dO' 



In 



h + ^/tlTi 



Ackno'wledginents 



L/2 y/E-COs20'' 

-fcBSn"^[sin(0 - 7r/2), ks] 

fc£;Sn~"'^[cOS0, fcs]. 



(65) 

(66) 

(67) 
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EXTENSION TO A PERIODIC AND FINITE 
WIDTH SOLUTION 

The analogy of the E-L equation [Eq. (fT8|) ] with a ficti- 
tious particle dynamics allows an extension of the single 
soliton solution to a periodic soliton array. The latter 
is obtained by choosing the integration constant E < 1, 
corresponding to the particle starting with a vanishing 
velocity and below the potential maximum. The subse- 
quent "evolution" of ^^(i) is clearly periodic in t, confined 
to the range < 9s{t) < 9m, with 9^. = |arccot(-E). 

Going back to Eg. ([TO]) we observe that the 9' integral 
can be related to the Legendre form of the elliptic integral 
of the first kind. 



r4 

F{<f>,k)^ / 
Jo 



^ 



(60) 



where fc^; — ■y/2/(l -t- E), and we used the fact that 
sn[(/), fc] is an odd function of </>. Thus the periodic conical 
soliton solution is given by 



9,{t,k) 



arccos < sn 



llnfi±^^ 
fc Vto + \/tg + l 



(68) 



with fc = 1+ giving our earlier single soliton solution 
[Eq.dUD]. 

One might hope to use this solution to model a fi- 
nite thickness, w, symmetric liquid- crystal cell with two 
boundaries inducing a symmetric (about z = u'/2) direc- 
tor rotation from 9 = 7r/2 to and back to ti/2. Naively, 
this maybe done by choosing the value of fc such that 
the period matches the cell thickness, w. Although this 
is possible for standard ID solitonic problems, because 
here the solution is periodic in t = [z + ZQ)/r (stemming 
from the fact that we are dealing with a 2D problem) and 
not in z, solution (|68p cannot be used to model a cell with 
sym,m,etric boundaries at fixed z = and z = w. A more 
general class of solutions is necessary but is currently un- 
available. 
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